arXiv:1503.05926v2 [astro-ph.HE] 31 Mar 2015 


Mon. Not. R. Astron. Soc. 000,^?? (2015) Printed 1 April 2015 (MN style file v2.2) 

Extremely narrow spectrum of GRB110920A: further 
evidence for localised, subphotospheric dissipation 

S. F. Ryde^’^, B. Ahlgren^’^, J. M. Burgess^’^, J. Larsson^’^, A. Pe’er^, 

C. Lundman^, M. Axelsson^’^’®, S. McGlynn® 

^Department of Physics, KTH Royal Institute of Technology, AlbaNova University Center, SE-106 91 Stockholm, Sweden 

^ The Oskar Klein Centre for Cosmoparticle Physics, Stockholm, Sweden 

^Department of Physics, Stockholm University, Sweden 

'^Erasmus Mundus Joint Doctorate in Relativistic Astrophysics 

^Department of Physics, Columbia University, New York, USA 

® University College Dublin, Ireland 

^ University College Cork, Ireland 

^Department of Astronomy, Stockholm University, Sweden 


Accepted... Received...; in original form ... 


ABSTRACT 

Much evidence points towards that the photosphere in the relativistic outflow in GRBs 
plays an important role in shaping the observed MeV spectrum. However, it is unclear 
whether the spectrum is fully produced by the photosphere or whether a substantial 
part of the spectrum is added by processes far above the photosphere. Here we make 
a detailed study of the 7 —ray emission from single pulse GRB110920A which has a 
spectrum that becomes extremely narrow towards the end of the burst. We show that 
the emission can be interpreted as Comptonisation of thermal photons by cold elec¬ 
trons in an unmagnetised outflow at an optical depth of r ^ 20. The electrons receive 
their energy by a local dissipation occurring close to the saturation radius. The main 
spectral component of GRB110920A and its evolution is thus, in this interpretation, 
fully explained by the emission from the photosphere including localised dissipation 
at high optical depths. 

Key words: gamma-ray burst: general — radiation mechanisms: non-thermal — 
radiation mechanisms: thermal 


1 INTRODUCTION 

Photospheric emission has currently become the pivot of 
the study of radiation mechanisms in gamma-ray bursts 
(GRBs). In the relativistic fireball model (e.g. Meszaros 


(20061), large amounts of energy is injected at the base of 
the flow where the optical depth is huge and the photon field 
gets efficiently thermalised. As the fireball expands the op¬ 
tical depth of the plasma eventually reaches unity and the 
photons are able to decouple from the plasma. The region 
from where the photons escape is called the photosphere. 

Two main perspectives of photospheric emission mod¬ 
els are currently discussed: In the first perspective, the ob¬ 
served spectrum is the result of emission from two different 
emission zones, the photosphere and an optically-thin region 


(Meszaros & Rees|2000 Zhang & Meszaros|2002 Daigne & 

■■ ■ i-- “ ■ ■ 


Mochkovitch 20021). The emission from the photosphere is 


assumed to be close to a blackbody, implying a passive flow 
below the photosphere. The dissipation of the kinetic en¬ 
ergy (or alternatively Poynting flux) that takes place in the 
optically-thin region results in non-thermal processes such as 


synchrotron emission (Tavani 19961 and/or inverse Compton 
scattering ( Ghisellini et al.|2000 ) producing the non-thermal 
part of the spectrum ( Mfezaros et al.||2002 Deng & Zhang 
2014). Based on such an interpretation, observed MeV spec¬ 


tra have been fitted with a blackbody and a power law (over 


the CGRO BATSE energy range; Ryde (2004 2005 I, Ryde fc 
Pe’er| ( |2b09 1) or a blackbody combined with a Band function 
(over the Fermi energy range Guiriec et al. (20111; Axelsson 
et al. (2012); lyyani et al. (2013); Guiriec et al. (2013)). 


In the second perspective, the entire spectrum is due 
to emission in the vicinity of the photosphere; additional 
processes cause the spectrum from the photosphere to differ 
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from a blackbody (Rees & Meszaros 2005). Energy dissi¬ 
pation below the photosphere (subphotospheric dissipation) 
can cause significant broadening of the thermal component 
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2 S. lyyani et al. 


in the flow. The dissipation can be imagined to be local, at 
a certain position in the flow (Rees fc Meszaros|2005 Pe’er 


et al. 2005 Beloborodov 2010[|, or be contin uous, having 
an effect through out the flow ||Giannios|2012 Beloborodov 


2013). Moreover, the photosphere is, in general, not ex¬ 


pected to be a surface, defined by a single photospheric 
radius, but rather an extended volume. Variations in the 
temperature due to expansion (Beloborodov 2011) as well 
as high latitude effects on the Lorentz boost of the emis¬ 
sion ( Abramowicz et al.|1991 1 will cause the observer to see 
a multi-temperature emission, significantly broadening the 
spectrum (Pe’er 2008, see also Goodmanl986, Lundman et 
al. 2013). 

First observational evidence for photospheric, spectral 
broadening was given by GRB 090902B, during which there 
was an onset of a broadening mechanism: Initially the spec¬ 
trum is well described by a (narrow) multicolour blackbody, 
not much broader than a blackbody (Ryde et al. ]Mol | |Zhang| 


et al. 2011). Later, the burst spectrum evolves and gets 


broader, now resembling a more typical (broad) spectral 
shape (Ryde et al. 2011). In this paper, we study the 160 
s long pulse in GRB110920A, which is very well suited for 
time resolved spectral analysis. Initial analyses ( |McGlynn| 
et al.|20l2 Shenoy et al.|20r3 1 have shown that the spectra 


are hard and part of the spectrum has been suggested to be 
of photospheric origin. Therefore, this burst is a strong can¬ 
didate to detect the photosphere and to study its properties. 
We find that the spectrum and its evolution gives additional, 
strong support for the existence of energy dissipation below 
the photosphere, shaping the observed spectrum. 

This paper is organised as follows: describes the ob¬ 

servations and results of the spectral analysis of the burst; 

discusses the physical scenario resulting in the observed 
spectra; ij^is a general discussion about the various proper¬ 
ties of the burst and finally gives the conclusion. 


2 OBSERVATIONS AND SPECTRAL 
ANALYSIS 

GRB110920A was observed on 20 September 2011 by 
the Gamma Ray Burst Monitor onboard Fermi Gamma- 
ray Space Telescope . It has a fluence of 1.74 ± 1.24 x 
10“‘^ergcm“^ in the energy range 10 keV - 40 MeV. There 
is no detection of Large Area Telescope (LAT, 100 MeV - 
300 GeV) and LAT Low Energy (LLE, 30 - 130 MeV) emis¬ 
sion for this burst. No afterglow has been detected for the 
burst and therefore, the redshift, z, is unknown. The light 
curve of the burst is a single pulse with a Tgo = 161 s, 
the time during which 90% of the emission is received. The 
duration is much longer than the average long burst, with 
only ~ 5% of bursts having a longer duration (von Kienlin 


et al. 2014). The variability timescale is much longer than 
the estimated dynamical time of the burst which is of the 
order of tdyn ~ 0.2 ms (see (2.2.3). Such smooth pulses en¬ 


able us to follow the spectral evolution more distinctly in 
comparison to other bursts where there are many spikes in 
the light curve which makes it difficult to resolve the differ¬ 
ent episodes of emission. Figure[^shows the composite light 
curve of the sodium iodide (Nal) and bismuth gallium oxide 
(BGO) detectors of GBM, that detected the burst (Meegan 



Figure 1. The light curve of GRB110920A as was observed by 
the GBM detectors: the brightest sodium iodide detector, Nal 3 
(8 - 800 keV; upper panel) and BGO 0 (200 keV - 40 MeV; lower 
panel) are shown. 


et al.|[2009 ). The energy flux in the energy range 10 keV - 
40 MeV peaks at 12 s with 5.8 ± 0.55 x 10“®ergcm“^ s“^. 

A standard time resolved spectral analysis of the 
burst with HEASARG’s spectral analysis software package 
XSPEG (Arnaud 1996), using PHA (Pulse Height Anal¬ 
yser) data from Nal 0,1,3 and the BGO 0 detectors in the 
energy range 10 keV - 40 MeV was performed. The burst 
was analysed for its duration of 0 to 200 s. After 60 s there 
was no significant BGO emission. However, the BGO was 
included in the analysis from 60 to 200 s in order to include 
the upper limits of BGO emission. The time bins for the 
spectral analysis were chosen for a signal-to-noise (S/N) ra¬ 
tio = 60, thereby enabling us to constrain the evolution of 
the spectral properties with high level of significance (see 
Burgess] ( (MIt ). 

In search of the spectral shape of the emission, three 
empirical functions for the photon flux {Ne{E), in units of 
photons/cm^/s/keV) are used to fit the spectrum: 

(i) The Band et al. (1994) function, which consists of two 
smoothly joined power-laws, and has four free parameters: 
Epeak, the peak in the uFi, spectrum, a the asymptotic 
power law index below the peak and /3 the asymptotic power- 
law index above the peak, apart from the normalisation. 

(ii) The blackbody (BB), which has two free parameters: the 
temperature, T, and the normalisation, A; 


Ne{E) = A 


explE/kT] - 1’ 


( 1 ) 


where k is the Boltzmann constant. 

(hi) The power-law (PL), which has two free parameters: 
the power-law index, s, and the normalisation, K: N{E) = 
K(E/Eoy, where Eq = 1 keV. 

The Band function and power law have no direct phys¬ 
ical meaning, however the Band function can be interpreted 
as being the result of an emission mechanism, such as syn¬ 
chrotron or inverse Compton emission. Likewise, the power- 
law can be a valid approximation of various emission pro¬ 
cesses over a limited energy range. The blackbody can be 


1 
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GRB110920A: extremely narrow spectrum 3 


interpreted as the thermal emission from the photosphere of 
the outflow. 

In the following, we use different combinations of these 
functions to asses three general models: a single non- 


thermal emission component (Band function alone, (2.1.11, 


a two emission-zone model (Band-|-BB, (2.1.2 I and finally a 
Comptonisation model, which is characterised by two tem¬ 
peratures, the injected photon temparature and the elec¬ 
tron temperature. The latter model also includes a power 
law component which can be associated to an optically-thin 
emission, not directly associated with the main Comptoni¬ 


sation event (BB-(-BB+PL, (2.2|. 


The high-energy photon-index, /3, is very hard in most bins 
and mainly has upper limits. 

Figure shows that the Band Speak decreases from 
nearly 700 keV to 50 keV, following a broken power law with 
a break at ~ 28±3 s. We note that this break does not coin¬ 
cide with the flux break at ~ 12 s. The temporal power-law 
index before (after) the break is —0.35±0.07 (—0.95±0.02). 


2.1.2 Band + blackbody fits 

Multicomponent spectra for GRBs were first suggested by 


Meszaros & Rees (20001, where, e.g., a blackbody models the 


2.1 Inability of the Band function to fit the data 

GRB spectra are typically fitted with standard empirical 
models consisting of a single Band function with or without 
an additional blackbody. For GRB110920A such functions 
do not produce the best fit to the data as will be shown 
below. 


2.1.1 Single Band function fits 


photosphere and a non-thermal emission is expected from an 
optically-thin emission process (two emission-zone model, 
see also M&zaros et al. ( 2002[ ); Zhang & Meszaros (20021; 
Daigne & Mochkovitch ( 2002| . Indeed, |Ryde ( 2004[ 20051 


successfully fitted such a model to data from the CGRO 
Burst and Transient Source Experiment (BATSE). Due to 
the limited energy range of the detector (~ 20 — 2000 keV), 
the model was limited to a blackbody and a power-law func¬ 
tion modelling the non-thermal emission. Similarly, Fermi 
observations (8 keV - > 40 MeV) have been analysed with 
such a model in mind, where a blackbody is instead corn- 


A Band function was fitted to the spectra of individual time 

bined with a Band function (e.g., Guiriec et al. (|2011l. Ax- 

bins through out the burst. This is the standard analysis 

elsson et al. (2012|), or synchrotron emission (|Burgess et al. 

that is performed in, for instance, the Fermi GBM (Gru- 


2014[|Yu et al.|2015|). The Band (or synchrotron) model adds 

ber et al. 20141 and CGRO BATSE catalogues (Goldstein 

a second spectral peak at higher energies, producing double 

et al.|2013 1. The fits are usually interpreted as a single non- 

humped spectra. 


thermal emission component, such as synchrotron or inverse 
Compton emission. 

A sample spectral fit for GRB110920A is shown in the 
upper left panel in Figure For an acceptable fit the resid¬ 
uals are expected to be random. It is thus obvious from the 
wavy structure of the residuals that the Band function alone 
does not reproduce the actual spectral distribution of the 
data. This fact strongly suggest that the spectrum is more 
complicated than a Band function alone. Additional spec¬ 
tral breaks can be imagined as well as additional spectral 
components. 

The main problem with these fits is the inability of the 
low-energy power-law to capture the additional curvature, 
while the peak energy is well determined. The Band function 
fits do, however, indicate that the spectra are exceptionally 
hard and narrow (see also ^2.3[ ). The determined photon in¬ 
dex, a, is found to lie above -0.4 through out the burst, which 
is inconsistent with both slow cooling synchrotron (q ^ - 
0.670 as well as fast cooling synchrotron (a ^ -1.5) emis¬ 
sion processes. We also find that a increases linearly with 
time from —0.4 to ±0.8. This is contrary to the typically ob¬ 


Band function and blackbody was fitted to the spectra of the 


individual time bins through out the burst (see also McG- 
|lynn et ar| ( |2012| ; |Shenoy et al.| ( |^13| ); |Burgess et aL] | |2014[ ). 
This model is a better fit with respect to the Band function 
model as the residuals become more random (upper right 
panel in Figure]^. Moreover, the fits significantly improve 
the pgsta0by > 25 in certain time bins. In Figure^we show 
the difference in pgstat, with respect to the Band function 
fits, for different models used for the analysis. 

The fit results (Band Epeak and kT) are shown in Fig¬ 
ure 0 When a blackbody is added to the Band function, 
the blackbody captures the low-energy spectral break and 
the Band Epeak gets pushed to higher values compared to 
the peak values found in the Band-only fits. This results in 
the residuals becoming more random which implies a bet¬ 
ter fit to the data. The presence of a blackbody leads to 
the low-energy photon-index, a, becoming softer, as is gen¬ 
erally observed (see, e.g., [Guiriec et aT](|2013|), but see also 


Burgess et al. (20141). During the burst, Epeak is observed to 


served behaviours where a evolves from hard to soft (Crider 


et al. 1997 Kaneko et al. 2006 Ghirlanda et al. 20031. In 


particular, after 70 s, a lies above 0 making it incompati¬ 
ble with any optically-thin, non-thermal processes and these 
spectra are among the hardest GRB spectra ever observed 


(Ryde et al. 2010 Kaneko et al. 20061. Moreover, we find 


that Epeak and a have a negative correlation: Epeak decays 
almost exponentially with respect to a, which is clearly op¬ 


posite to what is typically observed (Kaneko et al. 20061. 


^ The expected measured value for slow cooled synchrotron emis¬ 
sion is even softer, closer to a ~ —0.8 (Burgess et al. 2014) 


decrease from nearly 2 MeV down to approximately 20 keV. 
The Band function narrows with time and towards the end of 
the burst, the spectrum around the peak becomes nearly 
as narrow as a blackbody. The low-energy power-law index 
a increases from -1 to values larger than 1. It is worth noting 
that in this model also, after 100 s, a becomes harder than 
0. However, we note that before 100 s, the value is nearly 
constant at a ~ — 1 which is consistent with slow cooling 
and (modified) fast cooling synchrotron emission. The high- 
energy photon-index, /3, is very hard in most bins and has 
only upper limits. The correlation between the Band Epeak 
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Figure 2. Spectral fits to the time interval 11.28 s —15.2 s, green/solid line for Nal detectors and blue/solid line for BGO detector, (a) 
Band function fits. The residuals are wavy close to the peak of the spectrum and it clearly does not fit the data at lower energies, (b) 
Band + BB fit, the wavy structure of the residuals is less pronounced, but still remains, (c) The best fit model, Comptonised (2 BB) + 
power law. 


and a does not show a clear negative correlation, since a is 
nearly constant until 100 s. 


The blackbody temperature, kT, decreases smoothly 
with time with no significant break (Fig. [^. This is in 
contrast to what is typically observed where the tempera¬ 
ture decay has a characteristic temporal break (|Ryde|200^ 


Ryde fc Pe’er|2009| [Axelsson et al.|2012[ [Penacchioni et ah] 


2012[ Burgess et al. [2014 1. The relative strength of the non- 

thermal component decreases with time: ratio of the black- 
body flux to the total observed flux, Fbb /Ftot, evolves from 
~ 30% to nearly 90%, which is larger than typically observed 
values (Burgess et al.pOldb 


Alternative fits for a two-zone model for GRB110920A 


have been performed by Burgess et al. (20141; Yu et al. 


(20151, who replaced the Band function with an optically- 


thin synchrotron emission function. They find that only slow 
cooled synchrotron emission is allowed by the data, mainly 
due to the restrictive curvature of the spectrum around its 
peak. 


2.2 Best fit model: 2 blackbodies -|- power law 


As obvious from the above fits ((2.1.11, the spectra are more 
complex than a simple Band function. The shape of the spec¬ 
tra near the spectral peak is more of a top-hat; at the same 
time the spectra are very narrow. In addition to this, the 
low-energy break in the spectrum cannot be well captured 
by neither the Band function nor the BB -|- Band function 
fits. This motivates us to explore a Comptonisation model. 
Indeed, in Pe’er & Waxman ( 2004[ ) and Pe’er et al. 2006, 
it was shown that the process of Comptonisation can re¬ 
sult in such top-hat spectral shapes (see also discussion in 


Ghisellini & Celotti (1999l). 


2.2.1 Motivation of empirical model 

Assume that the photospheric emission, that is approx¬ 
imated by a blackbody with temperature Ti, undergoes 
Comptonisation by a thermalised pool of energetic electrons 
of temperature Tg. The resulting spectrum will then be char¬ 
acterised by two temperatures Ti and Te. Under certain cir- 
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GRB110920A: extremely narrow spectrum 5 


[h] 



Time (s) 


Figure 3. The Band function peak, -Fpeak? evolves as a broken 
power law from 700 keV to 50 keV with the break occurring at 
~ 28 s. Note that these fits are for the Band only model, which 
is not the best fit model. 
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the fitting routine becomes greatly simplified compared to 
a fitting routine invoking a fully developed Comptonisation 
model (see Ahlgren et al. 2015, in prep.) for the fits, but 
still retaining the important features. This is in particular 
important when exploring such spectral behaviour in large 
amounts of data. 

The two-blackbody approximation further neglects any 
broadening effects due to high-latitude emission (Pe’er 2008, 
Lundman et al. 2013). This is motivated by the fact that sev¬ 
eral observed burst spectra are inconsistent with the shape 
expected from the simplest models of such broadening. This 
is particular the case for the blackbody bursts in Ryde 


(2004), and more recently for GRB100507 Ghirlanda et al. 
(20131 and GRB101219B Larsson et al. ( 2015[ ). These bursts 
point towards the emitting region not being spherically sym¬ 
metric (which is typically the underlying assumption), but 
rather containing patches on scales smaller than l/T. Sim¬ 
ilarly small-patch emission zones are also found in the jet 
simulations (Lopez-Camara et al. 20131. In particular, as 


shown in §2.3[ the width for GRB110920A is also inconsis¬ 
tent with such broadening, in particular, towards the end of 
the burst. This argues that such effects should be secondary 
in GRB110920A. 

In §3.3| and 33.4| we illustrate how the two-blackbody 
approximation is used in interpreting the Gomptonisation 
process. 

With this motivation we use a fitting function that 
combines two blackbodies, that captures the main Gomp- 
tonised component, and a power-law function that, over 
the observed energy band, captures any secondary, broad¬ 
band emission. Such a power-law component was indeed 
found to accompany the photospheric emission component 
in GRB090902B (Ryde et al. 2010| and was suggested to 
be due to dissipation processes well above the photosphere 
(Pe’er & Ryde 20111. We note further that the fit func¬ 
tion consisting of two blackbodies and a power-law function 
has equal number of parameters as the Band + blackbody 
model, which simplifies model comparisons. 


Figure 4. The difference in pgstat for different models (Band-I- 
BB (black/squares) and Comptonisation (2BBs)-|- power law 
(red/circle)) with respect to the Band function fits done to the 
spectrum. Also, the difference in pgstat of the model Comptoni¬ 
sation (2BBs)-|- power law with respect to Band -|- BB model is 
shown (blue/triangle). 


cumstances the main features of such a spectrum can be 
captured by an approximation consisting of two blackbodies, 
with temperatures Ti and Te. We will show in lyyani et al. 
(2015 in prep.) that such an approximation is valid for ratios 
of Ts/Ti not much larger than ~ 10, optical depth less than a 
few tens and when there is no significant synchrotron photon 
production as a result of dissipation. Furthermore, an under¬ 
lying assumption of using the high-energy blackbody is that 
the electron distribution is Maxwellian. This is, in general, 
a good approximation for radiation mediated shocks below 
the photosphere (see, e.g., Bromberg et al. (20111). 

The advantage of using a simple analytical model is that 


2.2.2 Results of the empirical fits 

The model fits the data very well as can be seen in the 
example in Figure We find that until 100 s, this model 
is indeed a better fit with respect to the Band -|- BB fits 
and significantly improves the pgstat (by > 35 in certain 
time bins), see Figure]^ After 100 s, we find the power-law 
component is no longer statistically needed. Therefore, the 
spectra after 100 s are fitted only using two blackbodies, 
having merely four free parameters. 

In Figure we plot the evolution of the spectra of 
only the Gomptonised component (omitting the power-law 
component). In the figure the spectra are shifted in energy 
in order to highlight the evolution of the spectral shape. 

Since the model describes the shape of the spectrum 
well and also proves to be statistically significantly better 
than Band -|- BB fits in the beginning of the burst, we con¬ 
sider this model to be the best model for the burst. 
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6 S. lyyani et al. 




Figure 5. Results of Band function + blackbody fits: (a) The Band peak, Fpeak^ decreases from 2 MeV to 20 keV with no significant 
break, (b) The blackbody temperature, kT decreases smoothly and monotonously with time, in contradiction to the typical broken 
power-law observed. 



Figure 6. The evolution of the main spectral emission component of every second time bin is shown (Fu versus energy). Note that the 
time-resolved spectra are arbitrarily shifted in energy in order to highlight the change in spectral shape. 


2.2.3 Best fit parameters 

Properties of the fits are shown in Figurej^ The peak of the 
low-energy blackbody (BBl)Q J5i, decreases from nearly 200 
keV to 8 keV with a break at 12.9±1.7 s (Fig. Izt)- The break 
coincides with the peak in the light-curve. The power-law 
indices are before (after) the break —0.1±0.08(—0.69±0.06). 

This evolution is very different from the cooling be¬ 
haviour found in the Band-|-BB fits above, where the tem¬ 
perature did not show any break in its evolution. The break 
in the Ei evolution found here is, on the other hand, in agree¬ 
ment with what is typically observed for cooling blackbod- 
ies ( Ryde|2004 Ryde fc Pe’er|2009 Penacchioni et al.|2012 


Axelsson et al. 20121. This fact is indeed reassuring and 


strengthens the conclusion that the correct spectral model 
is used. 

The spectral peak of the high energy blackbody (BB2fl 
Ec, also decreases as a broken power-law from nearly 500 
keV to 40 keV, again with a break at 13.2 ± 2.4 s, which 
is consistent with the break in Ei. The power-law indices 
are before (after) the break —0.06 ± 0.19(—0.93 ± 0.03). We 
note that the post-break index is significantly steeper than 
for Ei = Ei{t). This is also illustrated by Figure]^, in which 
the ratio Ec/Ei is plotted. The model-ratio is shown as a 
blue line, has a maximum of ~ 3 and gets smaller with time. 


* Note that the thermal peak is given by E; = 2.7kTi where T; is 
the temperature of low-energy blackbody 


® Note that the thermal peak is given by Eq = 2.7kTg where Te 
is the temperature of high-energy blackbody 
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GRB110920A: extremely narrow spectrum 7 


This suggests a decreasing spacing between the peaks (see 
further the discussion in (4.21. Moreover, different from the 


evolution of Ec, after 50 s, Ei{t) does not show any clear 
trend, and does not follow the initial power-law decay, and 
exhibits a fluctuating nature. Thus, Ei is highly correlated 
with Ec only until 50 s. 

Apart from the main emission component, the power- 
law component also varies, with its index changing from -1.5 
to nearly -4, but it gets more poorly determined as the burst 
progresses. The power-law mainly accounts for the emission 
in the BGO as well as for the emission below 30 keV. How¬ 
ever, after 65 s, there is no more detected BGO emission 
and after 100 s, the emission below 30 keV becomes less sig- 
nihcant. The power-law component of the fits thus becomes 
insignihcant and the spectra are consistent with only the 
Gomptonised component (two blackbodies). 

The flux of the low-energy blackbody, Tbb, is around 
30% of the total observed flux, Ttot, while the flux of the 
high-energy blackbody , Ec, dominates the total observed 
flux through out the burst (Fig. IZli)- 

Finally, we note that the timescale of variations ob¬ 
served in the light-curve of the burst is much longer than 
the dynamical time which for the typical values of photo- 
spheric radius, rph = lO'^^ cm and Lorentz factor, F = 300, 
is given by rph/(2r^c) ~ 0.2 ms. This tells us that the cen¬ 
tral engine of the burst varies over a longer timescale and 
that the flow can be regarded as quasi-static in time. This 
validates the assumption that every time-bin can therefore 


out the entire duration of the burst, the SW is smaller than 
what is typically observed: the distribution of observed peak 
spectral widths has a maximum at SW ~ 10 (Axelsson et 
al. 2014, MNRAS, submitted). 

As an example, the Band function fit for the time in¬ 
terval 96 — 110 s is plotted in Figure together with a 
blackbody as well as the emission from a non-dissipative 
photosphere in a spherically symmetric outflow (broadened 
blackbody), aligned at the same uE^ peak energy. The spec¬ 
tral width of the Band function fit is SW = 5. Finally, we 
also studied the spectral width of the best fit model (Comp- 
tonisation model, (2.2|: The SW decreases with time from 


7 to 3.6 which is very similar to what was observed in the 
case of Band-only fits (Fig. &)• 


3 PHYSICAL SCENARIO 

3.1 Quasi-thermal Comptonisation 

As the best-fit spectra suggest, Comptonisation of a ther¬ 
mal photon component is a strong candidate in explaining 
the observed shapes. Quasi-thermal Comptonisation in vari¬ 
ous scenarios has previously been suggested to explain GRB 


be analysed assuming a steady flow. 

1999 

1. 



Rees & Meszaros (2005 I suggested that dissipation due 


2.3 Spectral Width 

Different emission mechanisms produce predictable widths 
of their emission spectra. The spectral width, SW, is here 
given by the ratio of the two extreme values of energy at 
which the vFv value of the spectrum is equal to half its 
maximum peak ( Axelsson fc Borgonovo|2015 1. For instance, 
a blackbody function is very narrow SW = 3.5, while the 
observed spectrum from a non-dissipative photosphere, as¬ 
suming a spherically symmetric, coasting wind has SW = 


to internal shocks below the photosphere can lead to Comp¬ 
tonisation of the thermal distribution of photons entrained 
in the outflow. Such dissipation can also be due to hadronic 
collision shocks ( Beloborodov|20TT| , or recollimation shocks 
occurring due to interaction between the jet and the star as 


7 (Lundman et al. 20131. This is since the spectrum is a 


multicolour blackbody, rather than a pure blackbody due to 
additional flux of photons emitted at high latitude^ In com¬ 
parison, synchrotron spectra are even broader. The spectral 
width of synchrotron emission from a mono-energetic elec¬ 
tron distribution is SW = 8.6. However, such a distribution 
is not generally expected. Fast cooling synchrotron emis¬ 
sion (power-law electron distribution) has SW ~ 370, while 
slow cooling synchrotron emission has SW = 145 (assuming 
the electrons distribution is a combination of Maxwellian -|- 


the jet traverses the envelope of the progenitor star (Lopez- 
Gamara et al.|[?013 Duffell fc MacFadyen|[20l4 1. In partic¬ 
ular, as the jet is ejected out of the star such shocks are 
expected ( Mizuta fc Ioka|2013 |. 

At large optical depths, r, due to the balance between 
the heating (direct Compton scattering) and cooling (inverse 
Compton scattering), the electrons attain a steady state and 


power law (Baring et al. 1995 Baring 2011 Summerlin & 
Baring|[2012 K 


We, therefore, measured the spectral widths from the 
flts. The spectral width of the Band-only fits decreases with 
time from nearly SW = 7 to 4. This illustrates that the spec¬ 
trum around the vE^, peak is very narrow; the SW even ap¬ 
proaches the narrowness of a blackbody function. Through- 


® High-latitude emission is here defined as the emission that orig¬ 
inates from angles larger than l/F to the line-of-sight of the ob- 


remain sub-relativistic, with values of 7/3 = 0.1 — 0.3 ( Pe’er 
et al. poo^ . The thermal photons that are entrained in the 
flow from the central engine will serve as seed photons that 
will be Compton scattered by the energetic electrons. Be¬ 
fore entering the dissipation site, these photons have been 
fully thermalised deep in the flow, either at the central en¬ 
gine itself or at a larger distance at which the flow starts to 
accelerate (the jet nozzle), denoted here by tq. No magnetic 
field is therefore needed to form the seed photons. In such 
a case, when the seed thermal photons scatter off the sub- 
relativistic electrons, the change in energy per scattering of 
the photons is small. At the same time, due to large num¬ 
ber of scatterings, rise, the photons diffuse from the thermal 
pool to higher energies and eventually form a new peak as 
a result of saturated Comptonisation. 

The spectrum that is formed through Comptonisation, 
at such a dissipation site, will be advected with the flow until 
it reaches the photosphere and the photons are emitted to 
the observer. This includes adiabatic cooling of the photon 
field. 
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Figure 7. Results of the best fit model (2BB + power law): (a) The peak energy of the low-energy peak, Ei, and the high-energy peak, 
Eq. The blue lines indicate the best fit to a broken power law function. Note the fit to Ei is only performed until 50 s, after which it 
starts to deviate, (b) The ratio EqJEw the solid line corresponds to the ratio of the best fit models in figure (a), (c) The spectral width 
of the vFn peak decreases with time from ~ 7 to 3.6, approaching the width of a blackbody. (d) The fraction of the observed total flux 
residing in the low-energy BB (FBe/Ftot) and in the high-energy BB (Fb/Ftot). 


3.2 Radial distribution of heating 

The spectral shape of the emission released at the photo¬ 
sphere in such a scenario will depend on the properties of 
the subphotospheric dissipation. Apart from the amount of 
kinetic energy that is converted, it will also depend on its 
properties such as the dissipation’s temporal variation and 
radial extent and its variation. For instance, even for steady 
flows, it can be envisioned that either the dissipation occurs 
at all radii with some prescribed radial dependence (con¬ 
tinuous dissipation) or that the dissipation is localised at 
a certain radial position (localised dissipation). Continuous 


dissipation can be due to radial and/or oblique shocks (Rees 


fc Meszaros|2005 Pe’er et al.|2006a Lazzati et al.|2009 l, or 


magnetic reconnection ( Giannios|2008| Begue & Pe’er p0l4| 

and in general leads to the formation of smooth Band-like 
spectra. All detailed spectral features produced by a dissipa¬ 
tion episode at a certain depth will, by necessity, be washed 
out by subsequent dissipations closer to the photosphere, 
where the jet properties have changed. Moreover, in the case 
of continuous dissipation, since efficient photon production 


is only expected deep down in the flow, the number-density 
of photons will not be able to maintain the level required 
for full thermalisation. The resulting peak (at high optical 
depths) will thus be a Wien peak with a temperature higher 
than a corresponding thermalised, blackbody spectrum. 

However, in a non-dissipative flow, even without pho¬ 
ton production, the spectrum remains thermalised and a 
blackbody distribution is retained. In such a case, a localised 
dissipation can distort the photon spectrum through Comp- 
tonisation of these blackbody photons and complex spectral 


shapes are easily achieved (e.g. Pe’er et al. (2006b|). Such 


localised dissipation can also be due to radial and/or oblique 
shocks, collisional processes | |Beloborodov|[MT^ , or due to 
some characteristic scale, e.g. the surface of the progenitor 
star ( Mizuta fc Ioka|2013 |. 

The existence of two distinct breaks in the spectrum of 
GRB110920A, and the fact that the spectra are extremely 
narrow, therefore strongly favours a scenario where the dis¬ 
sipation occurs over a range for which the jet parameters do 
not vary significantly, for instance, over a relatively narrow 
range in radii (localised dissipation). The dissipation site 


© 2015 RAS, MNRAS 000, [^?? 


















































GRB110920A: extremely narrow spectrum 9 



Figure 8 . Narrowness of the spectrum. The Band function fit 
for the time interval 96 s —110 s is shown in black (solid line). 
The blackbody (red/dash line) as well as the photospheric emis¬ 
sion from a spherically symmetric outflow at the coasting phase 
(blue/dash-dot line) is also plotted for the same peak energy to 
show the narrowness of the spectrum. The spectral width of the 
blackbody is SW = 3.5, the photospheric emission from a spheri¬ 
cally symmetric outflow is SW = 7 and that of the Band function 
in this time bin is 5. 


can therefore be characterised by an effective optical depth 
or, correspondingly, a dissipation radius, ru. Note, however, 
that any effect of the dissipation would occur over a range 
of radii determined by the timescales of the interaction be¬ 
tween the electrons and the photons. 

Furthermore, we argue that in the case of GRB110920A 
the dissipation radius, lies above, or close to, the satura¬ 
tion radius, rs. This is based on the fact that the dissipa¬ 
tion occurring in the flow alters the spectrum significantly. 
There should thus be significant kinetic energy in the outflow 
that could be dissipated and given to electrons. However, if 
rd < ra, then the kinetic energy density of the outflow would 
be less than the photon energy density, and a significant de¬ 
viation is not expected. Moreover, we do not expect to have 
dissipation during the acceleration phase. 


3.3 Seed photon component and outflow 
parameters 

We have argued that the flow in GRB110920A is a non- 
dissipative, passive jet, apart from a local, major dissipation 
episode at ra. The low-energy break in the observed spectra 
can thus be identified as the temperature of the seed photon 
distribution, that is advected from the central engine, and is 
injected in the dissipation site. This component allows us to 


calculate the outflow properties following Pe’er et al. (20071. 


However, many of the seed photons must have been 
scattered to higher energies due to the Gomptonisation pro¬ 
cess. Therefore the normalisation of the fitted (low energy) 
blackbody does not correspond to that of the injected seed 
blackbody. The photon number is conserved in the scattering 
process, therefore the seed photon normalisation can be esti¬ 
mated from the total number of photons in the Comptonised 
photon distribution. We note further that, additional pho¬ 


tons are not expected to be produced unless very deep In the 
flow ( Vurm et al.|2013 l or may be produced by synchrotron 
emission (Beloborodov 20131, which however depends on the 
strength of the magnetic field which is considered to be very 
low here (see (4.1 1 . Thus, the original number of photons 
in the fireball, Nq, is equal to the number of observed pho¬ 
tons, Nbbi + Nbb2, where Nbbi and Nbb2 are the number 
of photons in the two fitted blackbody components respec¬ 
tively. Nbb 2 corresponds to the number of photons scattered 
from the seed component. 

Using the estimated original number of photons we can 
calculate the Lorentz factor, F, the photosphere radius, rph, 
the saturation radius r^, and the nozzle radius of the flow. 
To. For these calculations to be valid we have to assume that 


the fra c tion o f kinetic energy dissipated is low (see Begue & 
lyyani (2014|) such that V = rj = Lq/Mc^^ where Lq is the 


total burst energy and M is the mass ejection rate, see fur¬ 
ther j ^3.4| Furthermore, the estimation of these parameters 
depend on the radiative efficiency of the jet {1/Y), which is 
unknown for GRB110920A. Here, Y is the ratio of the total 
burst energy to the observed prompt 7 — ray emission. How¬ 
ever, Racusin et al. (2011 1 estimated the radiative efficiency 


for a sample of 69 moderately bright bursts (with redshifts) 
and Y was estimated to be around 2 ^ Y 20 with an 
approximate average value of T ~ 5. In the following calcu¬ 
lations we therefore use T = 5. 

The estimated parameter values for a redshift, z = 2, 
are shown in Figure and [^. For instance, for the time 
bin at 10 s, the values of the outflow parameters are: 
F = 660±25, rph = 1.8±0.2x 10^^ cm , = 4.8±l.lx 10“ 

cm and ro = 7.3 ± 1.9 x 10^ cm. F decreases as a bro¬ 
ken power law, with a break at t = 15 ± 2 s with a 
slope before (after) the break given by the power law in¬ 
dex —0.04 ± 0.09 (—0.71 ± 0.04). The rph increases slightly 
until 50 s (rph oc after which it does not show 

any particular trend. While ro increases as ro oc ^0-923=0. 05 ^ 
rs exhibits a weaker increase (rs oc i0.4±o.04^^ note here 
that the estimated value of ro ~ 10® — 10® cm (Fig. [^) lies 
within the range of the predicted values for the formation 
of the first collimation shock that is found in numerical sim¬ 
ulations of hydrodynamical shocks moving though the en¬ 
velope of the progenitor star ( Duffell fc MacFadyen|[2014 1. 
The existence of such a shock is a stable prediction in these 
simulations, and the value of ro is typically around 10® cm, 
mainly depending on the dimensionless entropy of the flow. 
A high entropy flow (leading to high F) will have a smaller 
collimation shock radius and vice versa. 

The above calculations assume a classical fireball model 
(Pe’er et al. 2007). We neglect, for instance, effects of off-axis 
emission (high-latitude effects) which have been suggested 
to be an explanation of the observed behavior of the tem¬ 
perature and flux (Pe’er 2008[ Ryde & Pe’er 20091. Here, 
we instead assume that the varying observed properties re¬ 
flect varying properties of the flow at the central engine (di¬ 
mensionless entropy, temperature, and luminosity) and ra¬ 
dial (shock) dissipation pattern. For highly magnetised jets 
further alternatives exist by allowing the magnetisation at 
the central engine and magnetic dissipation pattern to vary 
( Hascoet et al.|2013 Gao fc Zhang|2(5l5 |. 

Finally, the outflow parameters are calculated for a red- 
shift, 2 = 2, an averaged value of GRBs, since the redshift 
of the burst is unknown. We also assume a fiat universe 
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Figure 9. (a) The evolution of the Lorentz factor, F, of the outflow, (b) The evolution of the photospheric radius, rph (black/squares), 
the nozzle radius, ro (blue/diamond), the saturation radius, Tb (green/triangles). The evolution of the radius (upper limit) of the localised 
dissipation, rj (red/circles). All the parameters are estimated for a T = 5 and a redshift, z = 2. 


= 0.73, i?o = 71). A diflerent value of z would only 
translate the values of the parameters depending on their 
dependences on a (for e.g. ^ = 0.3 lowers the estimated val¬ 
ues within a factor of 5) but will not change their temporal 
behaviour ( lyyani et al.|20l3 |. 

Summarising, the low energy component allowed us to 
determine the following flow parameters, F, Cph, rg, and ro. 


3.4 Comptonised peak and the dissipation site 


The Comptonised spectrum is expected to extend, in the 
comoving frame, from the energy of the seed photons £’) = 
2.7 kT' = 2.7 9' msC? up to the maximum energy that the 
photons can be up-scattered to £'^ = rrieC^/f, where 

7/3 is the electron momentum and the factor / takes into ac¬ 
count if the electrons are in the Thompson regime or not and 
has a value of 1 - ~ 3 { Pe’er et al.|2006b I; F/ and S'i are the 
temperatures of the Comptonised photons and of the seed 
thermal photons at the dissipation site, respectively. Here 
P = vjc, where v is the electron velocity, c is the velocity of 
light and 7 is the Lorentz factor of the electrons. 

For a relativistic Maxwell-Boltzmann distribution of the 
electrons, the observed spectral break at high energies, E^, 
therefore gives an estimate of the electron momentum. 


ilPf = 


E.f 

F TTleC^ 


2/3 


( 2 ) 


where the factor = (rph/ru)^^^ takes into account the 
adiabatic expansion from the dissipation site until the emis¬ 
sion is released at the photosphere. 

Furthermore, the critical Compton y-parameter is given 
by i/crit = \n{Ec/EP), where E^ and Ei are the observed 
electron temperature and thermal, seed temperature at the 
photosphere. Figure |10^ shows that i/crit remains close to 
1 through out the burst. The minimum number of scatter¬ 
ings, nsc.crit, needed to produce the spectrum can then be 
estimated as 


In/Ec/Si) 
4/3 (7/3)2 • 


( 3 ) 


Since t ~ rise for a relativistically expanding outflow, this 
corresponds to a lower limit to the optical depth Tmin = 
Usc.crit, which is plotted in Figu re [T 0 ] 3 . 

We can thus use equations ([^ and ^ to estimate Tmin 
and (7/3)min by making use of the estimate of F from the 
analysis of the seed photon component in (^3.3[). For in¬ 


stance, for the time bin at 10 s we find that Tn, 


20 Yr 


3/20 


and ('yP)'Lin ~ 0-03 at the dissipation sit^ During 

the whole burst duration, we find that the electron momen¬ 
tum, 7/3min is steady around ~ 0.17 and thus it is apparent 
that the electrons, at steady state, have subrelativistic ve¬ 
locities (cold electrons). The estimated value of Tmin corre¬ 
sponds to an upper limit of the dissipation radius ra.max = 
Tph/Tmin ~ 8 X 10^° cm. The dissipation thus occurs in 

the vicinity of the saturation radius (~ 5 x 10^° To"®cm) 


such that Tn 


<T < Tmax = (?^ph/rs) -- 


3/2 


The evolution of rd,max is included in Figure]^. Indeed, in 
order for Vd to lie above Tg the above calculations indicate 
that y ^ 4 giving a lower limit of the value of Y. 

The shape of the spectrum changes after 50 s, with the 
spectral peak becoming narrower, and suggests an increase 
in y > Uciit and correspondingly r > Tcrit, thereby increasing 
the strength of the BB2 component in the fits and creating 
the extreme, observed narrowness of the spectrum. The ex¬ 
act value of the j/-parameter must, however, be determined 
by detailed modelling of the spectral shape. However, we 
note that since the observed spectrum does have significant 
deviations from a blackbody (or a Wien spectrum), the dis¬ 
sipation should still be such that rj ^ Tg. Also, r is limited 
by Tph/rg oc Lq/T^. Moreover, r must be smaller than, say, 
100 in order to avoid the creation of a pure Wien peak. The 
observed spectral shape thus limits the possible values of r. 

An analytical expression for the electron momentum, 
( 7 / 3)2 dependencies on the flow parameters expected 

for subphotospheric dissipation in GRB flows is derived in 


Pe’er & Waxman (20051. Let us define Mph as the initial 


^ We use the notation = Y/10^ 
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energy density of the photon field and Mei as the injected 
energy density that the electrons have achieved due to the 
dissipation. Assuming that the main energy loss mechanism 
for the electrons is inverse Compton scattering that will bal¬ 
ance the heating due to energy dissipation of the flow kinetic 
energy, a steady state value for the electron momentum will 
arise 

Here the ratio of Mei and Uph can be estimated by us¬ 
ing the observed properties of the seed blackbody (i3.3|. 
If rs < rd < rph then, using the standard fireball theory 
( [Meszaros et al.|2002p and with Lq being the isotropic equiv¬ 
alent luminosity of the burst. 


- ^0 
dyrrlr^c 


(g) 


-8/3 


and 


Mel 


Lote 


dTrr^r^c 


— '^ph^e 


(g) 


2/3 


(5) 


( 6 ) 


where te is the fraction of the dissipated energy that goes 
into electrons, and noting that Vd = rp^lr, we have 


'^el 




(7) 


Combining equations 1^, Q, and 0 we can solve for 
the optical depth: 


9 rrieC^ 


3/7 


n3/7p-6/7 


(fcT) 




( 8 ) 


which can be estimated from the observables, Er, kT and the 


derived values of F, rph, and rs (from (3.31. This expression 


gives an analytical expression for the general dependencies 
and gives a rough estimate of the optical depth needed for 
a certain steady state electron temperature. As an exam- 


10 Ig 


15/28 


0.7 


pie the 10 s time bin in GRB110920A gives 
(ce = 0.1 was assumed) which is roughly consistent with the 
estimate from the width of the spectrum above. 

Similarly, noting that 

/rd\-2/3 

;) (g) <»' 


e' = 


/ Lo 


k 


rrisC^ KAnrsV^ca/ 


equations Q, Q, and 0 can be rearranged to yield the 
dependencies of the ratio 


Es 

E; 


1/2 - 
oc eJ T 


1/2 


(g) 


2/3 


^-l/8^y4rl/4 


( 10 ) 


where the dependence is strongest on r, Ce, and the ratio 
(rd/rs). Above we found that ra ~ rg and r to be moderate. 
Therefore, equation (101 implies that dissipation parameter. 


Ee, cannot be too large in order to explain the observed ratio 
of spectral peaks. We note further that, in the case that 
Td ~ Tg, equation 0 yields that Uei/Mph ~ te- 

The observed change in the shape of the spectrum after 
50 s, in which the high-energy peak becomes more promi¬ 
nent, requires the optical depths to increase, as discussed 
above. This will cause the steady state electron momentum 
(eq. 0 to decrease more rapidly than the observed temper¬ 
ature decrease. This in turn will cause the distance between 
the peaks to get smaller, as is indicated by the blue line in 


Figure m One possibility for the optical depth to increase 
is that, while Tph oc F“^ typically increases, one can imagine 
Td ~ 10^^ cm to be relatively steady (related to the progen¬ 
itor size scale). This will cause an increase in r = rph/rd. 

Summarising, the high-energy break allowed us to de¬ 
termine an additional parameter for the burst, namely the 
dissipation radius, rd. 


4 DISCUSSION 

Several different types of MeV spectra have been identi¬ 
fied in time resolved analysis of GRBs. A small fraction of 
the bursts are pure blackbody, others have a double peaked 
shape (often fitted by a Band -I- blackbody), while yet others 
are best fitted by a single Band function. 

One interpretation of the different types of spectra is 
therefore that, in addition to a blackbody (passive jet pho¬ 
tosphere), a substantial part of the spectrum is formed by 
processes far above the photosphere, e.g by optically-thin 
synchrotron emission. A varying combination of these com¬ 
ponents give rise to different spectral shapes. Another in¬ 
terpretation is that the spectrum is fully produced by the 
photosphere and the different shapes are due to varying dis¬ 
sipation properties below the photosphere (still allowing for 
dissipation above the photosphere, but not as a main ingre¬ 
dient). The analysis made above on GRB110920A supports 
the latter scenario. 


4.1 General scenario 


The top hat spectra observed in GRB110920A have a shape 
that is in between a pure blackbody spectrum, observed in 
CGRO BATSE | Ryde||2004 1 an d Fermi bursts IRyde et al. 
(20101; Ghirlanda et al. (20131; Larsson et al.f (20151 and 


the double-peaked spectra observed by Fermi, exemplified 
by GRB100724B (|Guiriec et al.|[MT| and GRB110721A 
(Axelsson et al. |20i2| |, where the ratio of the spectral break 
energies are typically EREi ~ 20 ( Burgess et al.|2014 l. The 
ratio of the spectral break energies in GRB110920A is, how¬ 
ever, much lower EREi ~ 3. 

Multiple spectral breaks in GRB spectra have been 
claimed by many studies. Using data from the PHEBUS 
experiment Barat et al. (1998) found that, apart from the 
typical spectral break at ~ 300 keV, an additional break ex¬ 
ists at around 1- 2 MeV. Similarly, |Ryde k, Pe’er (20091 ar¬ 
gues that the BATSE data (25-1800 keV) suggest additional 
peaks above the observed energy range (see also Battelino 
et al. pOOTt ), which was supported by observations by the 
EGRET instrument (Gonzalez et al. 20091. Moreover, Yu 


et al. (2014) A&A, submitted, found that a spectral model 
including two breaks fit the data equally well as the Band 
function in many cases. The ratio between the breaks was 
found to be typically around 10. The breaks Yu et al. (2014) 
introduced is for a synchrotron model, so the scenario is dif¬ 
ferent from what we argue for here. However, their study 
indicates that the single peak identified in Band function 
fits can be interpreted as covering several breaks. This thus 
suggests that multiple breaks in GRB spectra could be more 
common than previously thought. 

Furthermore, the Comptonised spectrum suggested 
above for GRB110920A, in the fits performed in this paper. 
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Figure 10. (a) The evolution of the critical Compton y parameter, ycrit- (b) The evolution of the lower limit of the optical depth at 
which dissipation occurs, ~ ^rit for a, Y = 5. 


is captured by an empirical model consisting of two black- 
bodies and a power law. Fits using a similar model have 
been performed on several bursts by|Basak fc Rao] (|2014[); 


Raghurama Rao et al. (20141, however, motivated by a dif¬ 


ferent scenario. In any case, these fits could be an indica¬ 
tion that a Compotonised spectrum indeed is common and 
that the two spectral breaks can be identified. |Basak fc Rao| 
(20141; Raghurama Rao et al. (20141 find the ratio between 
the spectral peaks to vary between ~ 3 and ~ 10 in their 
sample. 


Based on these observed spectral properties, it can 
therefore be speculated that all such spectra (varying ra¬ 
tios of the spectral breaks) are formed primarily due to a 
Comptonised spectrum from a localised dissipation. Further, 
secondary, effects on the detailed shape of the spectra can 
be caused by synchrotron photons produced at low optical 
depths, if there are strong enough magnetic fields (|Vurm| 
et al.|2013| ). 

The ratio between the breaks in the Comptonised spec¬ 
trum is given by equation ( |10[ ), with the strongest depen¬ 
dences on T and the ratio ra/rg. Since the spectral shape 
limits the possible variation in r, the distance between the 
spectral peaks are mainly set by ra/rg. The closer the dissi¬ 
pation radius is to the saturation radius the higher the tem¬ 
perature of the seed photons (low-energy break) will be since 
they are less affected by adiabatic losses. Similarly, the ratio 
of energy density in electrons to that in the photons Mei/uph 
will become lower due to lower adiabatic losses. This, in 
turn, will cause the electron momentum (high-energy break) 
to decrease (eq. |^. The distance between the two peaks will 
thus get smaller as the dissipation site approaches the satu¬ 
ration radius. 


Bursts with larger ratio between the peaks Ec/Ei, can 
thus be partly due to them having a different dissipation pro¬ 
file with a dissipation radius occurring substantially further 
away from the saturation radius. Indeed, for many bursts the 
ratio of derived rph/rg is typically much larger than what is 
observed in GRB110920A and have values rph/rg ~ 100 (e.g. 


Axelsson et al. (20121; Guiriec et al. ( 2011[); Burgess et al. 


(20141). This yields, in turn, larger values of ra/rg (assuming 


a certain value of r). Similarly, internal shocks are expected 
to occur with ra/rg = 2r ~ 300. Using such a value the 
Ec/Ei ratio increases substantially (eq. |10[|. 


4.2 Cooling of the seed photon distribution 

We pointed out above that the temperature of the seed pho¬ 
tons did not follow the expected power law after ~ 50 s (see 
Fig.l^r). It was further argued that, since the photon number 
is conserved in the Comptonisation process, the seed pho¬ 
ton distribution is smaller than it was before the dissipation 
episode. If the number of scattered photons is small the fit¬ 
ted temperature should not differ greatly from the original 
value, only the normalisation will decrease. However, if the 
number of scattered photons is large then the original dis¬ 
tribution will be lost and the fitted value of the temperature 
will not correspond to the original value. 

In Figure the photon flux in the different compo¬ 
nents is shown. Before 50 s the number of photons in the 
two blackbody components are similar. However, after 50 s, 
the Comptonised peak dominates by a factor of four. This 
thus supports the argument that the erratic behaviour of 
the seed photon distribution, after 50 seconds (as compared 
to the predicted smooth power-law decay) is due to strong 
Comptonisation (larger ^-parameter) leading to a loss of the 
possibility to measure the original temperature. It is inter¬ 
esting to note that for the second to last data point the 
fluxes are similar again and the measured temperature lines 
up with the power law decay, that was determined from the 
data points before 50 s. 


5 CONCLUSION 

We thoroughly studied the single pulse GRB110920A since 
it was a strong candidate for revealing the behaviour 
of the photosphere. Indeed, the extreme narrowness of 
the spectrum observed directly excludes non-thermal emis¬ 
sion models. We have demonstrated that the spectrum in 


© 2015 RAS, MNRAS 000, [^?? 












































GRB110920A: extremely narrow spectrum 13 



Figure 11. Photonfiux light curves for the total flux (black/ 
squares), the black body components (blue/ diamond and red/ 
circles) and for the power law component (green/ triangles). Pho¬ 
tons from the Comptonised peak dominated the emission after ~ 
50 s. 


GRB110920A is more complex than a Band function, hav¬ 
ing a top-hat shape with two spectral breaks. Such an obser¬ 
vation restricts the amount of subphotopsheric dissipation 
allowed as well as the radial extent of such disturbances. If 
there is a large range of radii where there is strong dissi¬ 
pation the spectral features such as breaks will be smeared 
out and the resulting spectrum would be a smooth Band 
like spectrum. We find a best fit model which includes lo¬ 
cal Comptonisation which is characterised by two tempera¬ 
tures, the seed photon temperature and the temperature of 
the electrons, in addition to a weaker power-law component 
(which is, though, only present during the first 100 s). The 
main Comptonised component enables us to determine the 
properties of the jet such as the bulk Lorentz factor, the 
photospheric radius and the dissipation radius. In particu¬ 
lar, we find that the optical depth of the dissipation site is 
at r ~ 20. The narrowness of the spectrum is mainly due to 
the dissipation occurring close to the saturation radius and 
the optical depth increasing over the burst duration. The 
main spectral component of GRB110920A and its evolution 
is thus fully explained by emission from the photosphere 
including localised dissipation at high optical depths. 

We argue that the top hat spectra observed in 
GRB110920A could be the missing link between the very 
narrow spectra observed in many bursts (in some cases 
even fitted by a pure blackbody, e.g. the Fermi bursts 
GRB100507 and GRB101219B) and the double-peaked 
shape observed in Fermi bursts, such as GRB100724B and 
GRB110721A.Thus, in the scenario of localised subphoto- 
spheric dissipation, the main distinction between various 
spectral shapes is based on where the localised dissipation 
occurs relative to where the flow saturates. More spectral 
fitting using a physical model including subphotospheric 
dissipation is needed to assess such a framework properly. 


ACKNOWLEDGEMENTS 

We thank the referee for useful comments on the manuscript. 
We acknowledge support from the Swedish National Space 
Board and Goran Gustafsson Foundation. SI is supported 
by the Erasmus Mundus Joint Doctorate Program by Grant 
Number 2011-1640 from the EACEA of the European Com¬ 
mission. 


REFERENCES 

Abramowicz M. A., Novikov I. D., Paczynski B., 1991, ApJ, 
369, 175 

Arnaud K. A., 1996, in Jacoby G. H., Barnes J., eds. Astro¬ 
nomical Data Analysis Software and Systems V Vol. 101 
of Astronomical Society of the Pacific Conference Series, 
XSPEC: The First Ten Years, p. 17 
Axelsson M., Baldini L., et al. 2012, ApJ, 757, L31 
Axelsson M., Borgonovo L., 2015, MNRAS, 447, 3150 
Baring M. G., 2011, Advances in Space Research, 47, 1427 
Baring M. G., Ellison D. C., Jones E. C., 1995, Advances 
in Space Research, 15, 397 
Basak R., Rao A. R., 2014, MNRAS, 442, 419 
Battelino M., Ryde F., Omodei N., Band D. L., 2007, 
in Ritz S., Michelson P., Meegan C. A., eds. The First 
GLAST Symposium Vol. 921 of American Institute of 
Physics Gonference Series, Simulation of prompt emission 
from GRBs with a photospheric component and its de¬ 
tectability by GLAST. pp 478-479 
Begue D., lyyani S., 2014, ApJ, 792, 42 
Begue D., Pe’er A., 2014, ArXiv e-prints, 1410.2730 
Beloborodov A. M., 2010, MNRAS, 407, 1033 
Beloborodov A. M., 2011, ApJ, 737, 68 
Beloborodov A. M., 2013, ApJ, 764, 157 
Bromberg O., Mikolitzky Z., Levinson A., 2011, ApJ, 733, 
85 

Burgess J. M., 2014, MNRAS, 445, 2589 
Burgess J. M., Preece R. D., Connaughton V., et al. 2014, 
ApJ, 784, 17 

Crider A., Liang E. P., Smith 1. A., Preece R. D., Briggs 
M. S., Pendleton G. N., Paciesas W. S., Band D. L., Mat- 
teson J. L., 1997, ApJL, 479, L39 
Daigne F., Mochkovitch R., 2002, MNRAS, 336, 1271 
Deng W., Zhang B., 2014, ApJ, 785, 112 
Duffell P. C., MacFadyen A. L, 2014, ArXiv e-prints, 
1407.8250 

Gao H., Zhang B., 2015, ApJ, 801, 103 
Ghirlanda G., Celotti A., Ghisellini G., 2003, A&A, 406, 
879 

Ghirlanda G., Ghisellini G., et al. 2013, MNRAS, 428, 1410 
Ghisellini G., Celotti A., 1999, ApJ, 511, L93 
Ghisellini G., Lazzati D., et al. 2000, MNRAS, 316, L45 
Gianni os D., 2008, A&A, 480, 305 
Gianni os D., 2012, MNRAS, 422, 3092 
Goldstein A., Preece R. D., Mallozzi R. S., Briggs M. S., 
Fishman G. J., Kouveliotou C., Paciesas W. S., Burgess 
J. M., 2013, ApJS, 208, 21 

Gonzalez M. M., Garrillo-Barragan M., Dingus B. L., 
Kaneko Y., Preece R. D., Briggs M. S., 2009, ApJ, 696, 
2155 

Gruber D., Goldstein A., Weller von Ahlefeld V., et al. 
2014, ApJS, 211, 12 


© 2015 RAS, MNRAS 000, [^?? 









14 S. lyyani et al. 


Guiriec S., Connaughton V., et al. 2011, ApJL, 727, L33 
Guiriec S., Daigne F., Hascoet R., Vianello G., Ryde F., 
Mochkovitch R., Kouveliotou C., Xiong S., Bhat P. N., 
Foley S., Gruber D., Burgess J. M., McGlynn S., McEnery 
J., Gehrels N., 2013, ApJ, 770, 32 
Hascoet R., Daigne F., Mochkovitch R., 2013, A&A, 551, 
A124 

lyyani S., Ryde F., Axelsson M., Burgess J. M., et al. 2013, 
MNRAS, 433, 2739 

Kaneko Y., Preece R. D., et al. 2006, ApJ, 166, 298 
Larsson J., Racusin J. L., Burgess J. M., 2015, ApJL, 800, 
L34 

Lazzati D., Morsony B. J., Begelman M. G., 2009, ApJL, 
700, L47 

Liang E., Kusunose M., Smith 1. A., Crider A., 1997, ApJL, 
479, L35 

Liang E. P., 1997, ApJL, 491, L15 

Lopez-Camara D., Morsony B. J., Begelman M. C., et al. 

2013, ApJ, 767, 19 

Lundman C., Pe’er A., Ryde F., 2013, MNRAS, 428, 2430 
McGlynn S., Fermi GBM Collaboration et al. 2012, in - 
Ray Bursts 2012 Conference (GRB 2012) Photospheric 
Emission in Fermi Gamma Ray Bursts 
Meegan C., Lichti G., Bhat P. N., et al. 2009, ApJ, 702, 
791 

Meszaros P., 2006, Reports on Progress in Physics, 69, 2259 
Meszaros P., Ramirez-Ruiz E., et al. 2002, ApJ, 578, 812 
Meszaros P., Rees M. J., 2000, ApJ, 530, 292 
Mizuta A., loka K., 2013, ApJ, 777, 162 
Pe’er A., 2008, ApJ, 682, 463 

Pe’er A., Meszaros P., Rees M. J., 2005, ApJ, 635, 476 
Pe’er A., Meszaros P., Rees M. J., 2006a, ApJ, 642, 995 
Pe’er A., Meszaros P., Rees M. J., 2006b, ApJ, 642, 995 
Pe’er A., Ryde F., 2011, ApJ, 732, 49 
Pe’er A., Ryde F., et al. 2007, ApJL, 664, LI 
Pe’er A., Waxman E., 2004, ApJ, 613, 448 
Pe’er A., Waxman E., 2005, ApJ, 628, 857 
Penacchioni A. V., Rufiini R., Izzo L., Muccino M., Bianco 
C. L., Caito L., Patricelli B., Amati L., 2012, A&A, 538, 
A58 

Racusin J. L., Oates S. R., Schady P., Burrows D. N., de 
Pasquale M., Donato D., Gehrels N., Koch S., McEnery J., 
Piran T., Roming P., Sakamoto T., Swenson C., Troja E., 
Vasileiou V., Virgili F., Wanderman D., Zhang B., 2011, 
ApJ, 738, 138 

Raghurama Rao A., Basak R., Bhattacharya J., Chandra 
S., Maheshwari N., Choudhury M., Misra R., 2014, Re¬ 
search in Astronomy and Astrophysics, 14, 35 
Rees M. J., Meszaros P., 2005, ApJ, 628, 847 
Ryde F., 2004, ApJ, 614, 827 
Ryde F., 2005, ApJ, 625, L95 
Ryde F., Axelsson M., et al. 2010, ApJL, 709, L172 
Ryde F., Pe’er A., 2009, ApJ, 702, 1211 
Ryde F., Pe’er A., et al. 2011, MNRAS, 415, 3693 
Shenoy A., Sonbas E., Dermer C., et al. 2013, ApJ, 778, 3 
Summerlin E. J., Baring M. G., 2012, ApJ, 745, 63 
Tavani M., 1996, ApJ, 466, 768 

von Kienlin A., Meegan C. A., Paciesas W. S., et al. 2014, 
VizieR Online Data Catalog, 221, 10013 
Vurm L, Lyubarsky Y., Piran T., 2013, ApJ, 764, 143 
Yu H.-F., Greiner J., van Eerten H., et al. 2015, A&A, 573, 
A81 


Zhang B., Meszaros P., 2002, ApJ, 581, 1236 
Zhang B.-B., Zhang B., et al. 2011, ApJ, 730, 141 


© 2015 RAS, MNRAS 000, [^?? 


